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LONG-TERM  GOALS 

The  US  Navy  has  been  developing  its  Multi-Static  Active  Coherent  (MAC)  acoustic  search  system. 
MAC  is  an  air  launched  acoustic  search  system  that  combines  a  newly  developed  coherent  source 
sonobuoy  with  a  field  of  receiver  sonobuoys.  With  the  introduction  of  the  coherent  source,  operational 
search  data  collected  during  MAC  exercises  are  well  suited  for  post-processing  to  infer  waveguide 
properties.  The  long  term  goal  of  this  research  is  to  investigate  the  feasibility  of  estimating  sedient 
acoustic  properties  using  data  collected  during  such  exercises. 

OBJECTIVES 

Modal  dispersion  data  are  used  in  inversion  schemes  that  estimate  sediment  properties  in  range 
independent  and  range-dependent  environments.  The  experimental  scenario  for  the  use  of  air  launched 
sonobuoys  to  extract  information  about  the  sediment  properties  in  the  area  covered  by  the  sonobuoys  is 
to  have  one  of  the  sonobuoys  broadcast  a  broadband  signal  with  a  large  set  of  buoys  that  receive  these 
signals.  These  receiving  sonobuoys  are  dispersed  over  a  wide  area.  In  such  an  experimental  set  up  the 
data  acquired  by  the  sonobuoys  is  used  to  estimate  the  sediment  properties  in  the  area  covered  by  the 
sonobuoys. 

APPROACH 

The  objective  of  the  current  work  is  to  investigate  the  feasibility  of  estimating  sediment  properties 
from  acoustic  data  acquired  by  a  set  of  sonobuoys  distributed  over  a  wide  area.  This  is  done  using  a 
realistic  environment  and  analysis  methods  for  estimating  the  sediment  properties  from  mode 
dispersion  data.  We  create  the  experimental  set  up  by  locating  a  number  of  receiver  sonobouys 
uniformly  distributed  over  an  area  and  a  sonobuoy  source  broadcasting  a  broadband  signal.  The  signals 
acquired  by  the  receiver  buoys  are  computed  and  analysis  of  this  signal  is  done  to  obtain  mode 
dispersion  data  and  subsequently  the  sediment  properties  of  the  region  between  the  source  position  and 
that  of  the  receiver  sonobuoy.  An  additional  constraint  is  the  frequency  of  the  source.  If  the  source 
frequency  is  low  then  the  size  of  sources  that  generate  the  signal  are  typically  large  and  hence  difficult 
to  be  launched  by  a  sonobuoy.  As  the  frequency,  increases  the  size  becomes  smaller  and  therefore  the 
source  frequency  should  be  high  with  adequate  bandwidth.  The  inversion  procedure  therefore  is 
attempted  over  a  range  of  frequencies. 
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WORK  COMPLETED 


A.  Simulation  of  environment 

The  estimation  of  sediment  properties  in  shallow  water  areas  is  difficult  because  of  the  bathymetry  in 
these  areas  are  complex  giving  rise  to  varying  water  column  depths  as  well  as  varying  sediment 
characteristics.  In  order  to  represent  such  an  environment,  an  area  off  the  coast  of  New  Jersey  for 
which  the  bathymetry  is  known  was  used  in  this  simulation.  Figure  1  show  the  region  that  was  used  in 
creating  the  data  for  this  simulation  together  with  the  locations  of  25  sonobuoys  and  the  source.  The 
locations  of  the  sonobuoys  were  obtained  assuming  that  they  were  uniformly  distributed  over  the 
defined  limits  of  the  region  as  indicated  by  their  Longitudes  and  Latitudes.  The  area  covered  by  the 
sonobuoys  was  then  divided  in  to  six  regions  as  shown  in  Figure  1. 


Division  of  the  region  in  to  6  sub-regions 


Longitude  (deg  W) 


Figure  1:  The  panel  on  the  left  shows  the  location  of  the  distributed  network  of  sonobuoys  and  the 
location  of  the  source.  The  figure  also  shows  bathymetry  of  the  area  which  varies  from  60m  to 
100m.  The  right  panel  shows  the  six  major  regions  in  to  which  the  area  has  been  sub-divided. 

Because  of  the  differences  in  the  water  column  depths  along  the  different  source/buoy  transects  we 
made  assumptions  that  the  sediment  properties  also  varied  along  these  different  sections.  In  regions 
where  there  was  a  large  variation  in  water  column  depths  we  assumed  the  the  bottom  properties  were 
different.  For  example  in  the  area  where  Buoys  2,  3,  19,  and  20  are  located  there  are  three  distinct 
depth  regions  namely  region  of  an  average  depth  of  72  m  followed  by  a  region  of  an  average  depth  of 
69  m  and  then  the  final  section  with  an  average  depth  of  60  m.  In  this  case  we  assumed  that  the  regions 
where  the  depth  are  72  m  and  69  m  the  bottom  properties  are  identical  where  as  the  region  where  the 
depth  changed  to  60  m  has  a  different  bottom  properties.  Similarly  by  looking  at  the  depth  variations 
in  the  different  parts  we  created  a  total  of  seven  regions  with  slightly  different  sediment  properties.  In 
each  of  these  sections  different  bottom  models  were  assumed  and  these  are  termed  model  Ml  to  M7  as 
shown  in  Figure  2.  The  number  of  sediment  layers  and  their  thicknesses  have  also  been  varied  and  are 
similar  to  sediment  layering  found  in  the  area  in  earlier  experiments.  The  details  of  the  compressional 
wave  speed  profiles  for  the  seven  bottom  models  are  shown  in  Figure  2.  In  all  cases  the  terminating 
half  space  is  assumed  to  have  a  compressional  wave  speed  of  1850  m/s.  The  sediment  density  values 
were  set  at  1.6  gm/cc  in  all  the  layers  and  in  the  terminating  half  space.  Attenuation  in  the  sediment 
layers  and  terminating  halfspace  have  been  neglected.  The  sound  speed  structure  in  the  water  column 
was  assumed  to  be  the  same  over  the  entire  area. 
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Degtj^amd  bottom  models  for  the  different  parts  of  the  six  regions 


Compressional  Wave  Speed  (m/s) 


Compressional  Wave  Speed  (rtVs) 


Figure  2:  The  left  panel  of  the  figure  shows  the  six  major  divisions  of  the  region  together  with  the 
variations  in  the  depths  within  these  regions.  It  also  indicates  the  seven  bottom  models  used  to 
describe  the  range  dependent  sediment  properties.  The  compressional  wave  speed  profiles  of  the 
seven  bottom  models  are  shown  in  the  right  panel  of  the  figure. 

B.  Broad  Band  Signal  Generation 

Broadband  signals  that  have  found  wide  acceptance  are  the  LFM  (Linear  Frequency  Modulated) 
signals  or  chirps  in  which  the  frequency  is  varied  linearly  either  going  up  or  down  with  the  signal 
duration.  Broadband  signals  are  also  created  using  pseudo  random  sequences  [1],  The  sequence  of 
zeros  and  ones  generated  by  the  sequence  is  used  to  modulate  a  carrier  signal.  For  creating  the 
transmitted  signal  the  zeros  in  the  sequence  is  changed  to  -1,  (i.e.  we  now  have  a  sequence  of  1  and  -1). 
Each  digit  in  the  sequence  is  a  chip.  If  the  sequence  has  N  digits,  the  signal  has  N  chips.  The  length  of 
each  chip  defines  the  bandwidth  of  the  signal.  The  numbers  of  cycles  of  the  carrier  signal  that  form  a 
chip  determines  the  bandwidth  of  the  signal. 

C.  Computation  of  Signals  Acquired  by  the  Sonobuoys. 

In  order  to  compute  the  signal  acquired  by  the  buoys  the  distances  from  the  source  to  the  receiver  are 
required.  These  are  obtained  from  the  known  positions  of  the  source  and  the  respective  buoys.  A  time- 
domain  acoustic  signal  at  the  receiver  is  obtained  by  Fourier  synthesis  of  the  frequency  domain 
solution  to  the  acoustic  wave  equation  using  the  expression  below. 

p(r,z,t)  =  J^5((X»)p(r,z,  o))  exp(ia)t)  do)  (1) 

where  5  (m)  is  the  spectrum  of  the  transmitted  signal  and  p(r,  z,  m)  is  the  frequency  domain  solution  at 
range  r  and  depth  z.  To  compute  the  time  domain  solution  above  requires  that  p(r,  z;  m)  be  calculated 
at  a  number  of  discrete  frequencies,  which  span  the  frequency  band  of  the  source  spectrum.  The  ocean 
environment  was  the  input  to  the  normal  mode  code  KRAKEN  [2]  for  computing  the  frequency 
domain  solution  with  the  range  to  the  receiver,  depths  of  the  source,  and  receiver  set  as  required. 

The  signal  acquired  by  the  receiver  is  a  convolution  of  the  impulse  response  of  the  channel  with  the 
transmitted  pulse.  Before  processing  the  signal  acquired  by  the  buoys,  it  is  necessary  to  extract  the 
impulse  response  of  the  channel,  to  take  out  the  time-frequency  variation  of  the  transmitted  signal 
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during  signal  duration.  This  is  done  by  de-convolution  of  the  received  signal.  In  the  simulation  study 
presented  here,  the  deconvolution  of  the  signal  is  done  in  the  frequency  domain.  The  impulse  response 
is  then  used  to  perform  the  time-frequency  analysis  using  short  time  Fourier  transform  (STFT)  or 
warping  procedure  [3]. 

D.  Inversion  for  Compressional  Wave  Speed  Profiles 

The  inversion  for  the  sediment  compressionalwave  speed  profdes  for  the  seven  areas  as  represented  in 
models  Ml  to  M7  was  performed  using  the  mode  dispersion  data  obtained  from  the  signals  recived  at 
the  buoys.  The  procedure  for  carrying  out  the  inversion  is  well  documented  in  literature  [4,  5].  All 
inversions  were  carried  out  for  the  compressional  wave  speed  in  the  layers.  All  other  parameters 
namely  the  compressional  wave  speed  in  the  terminating  half  space  and  the  density  value  of  the 
sediment  layers  and  the  terminating  half  space  were  assumed  to  be  known. 

Consider  the  bottom  model  M2.  In  this  case,  the  water  column  has  a  depth  of  60  m.  The  sediment 
consists  of  three  layers  with  layer  thicknesses  of  8  m,16  m,  and  8  m  respectively.  The  compressional 
wave  speed  in  the  layers  are  1530  m/s,  1660  m/s,  and  1490  m/s  respectively.  The  terminating  half 
space  has  value  of  1850  m/s.  The  density  in  the  sediment  layers  and  the  terminating  half  space  is  1.6 
gm/cc.  Examination  of  the  mode  function  indicates  that  in  order  to  estimate  the  bottom  model  correctly 
we  need  to  know  the  mode  travel  time  at  frequencies  including  those  that  are  below  the  Airy  phase.  In 
many  cases,  the  strength  of  the  modes  at  frequencies  below  the  Airy  phase  is  very  low  and  therefor  it  is 
difficult  to  estimate  the  mode  travel  time  at  these  frequencies. 

In  order  to  estimate  the  values  of  the  mode  travel  time  at  frequencies  below  the  Airy  phase  the  bottom 
model  was  first  estimated  using  the  mode  dispersion  data  at  frequencies  above  Airy  phase.  Using  the 
bottom  model  thus  obtained  mode  dispersion  was  computed  theoretically  to  get  approximate  mode 
travel  times  at  frequencies  below  the  Airy  phase.  This  was  then  incorporated  in  the  input  data.  Using 
this  approach  repeatedly,  the  estimates  of  mode  arrivals  at  frequencies  below  the  Airy  phase  were 
incrementally  improved  and  the  more  correct  estimates  of  the  bottom  compressional  wave  speed 
profiles  were  obtained. 

Another  issue  that  needs  attention  is  determination  of  the  number  of  sediment  layers  and  their 
thickness.  If  the  number  of  layers  and  their  thicknesses  are  known  then  they  are  input  to  the  inversion 
process.  The  inversion  involves  obtaining  the  inverse  of  a  matrix.  This  matrix  is  generally  ill 
conditioned  and  need  some  form  of  regularization.  We  have  used  qualitative  regularization  [6,  7]  in 
obtaining  the  inverse  solution  and  this  requires  data  on  sediment  layering  as  input  to  the  inversion.  In 
order  to  have  an  estimate  of  the  number  of  layers  and  their  thicknesses  we  performed  inversion  without 
any  layer  data  (i.e.  without  incorporating  qualitative  regularization).  The  result  of  inversion  in  the 
region  of  model  M2  is  shown  in  Figure  3.  This  indicates  that  the  compressional  wave  speed  in  the 
sediment  has  a  low  speed  top  layer  which  is  followed  by  a  high  speed  layer  and  a  low  speed  layer 
again.  The  high  speed  layer  in  the  middle  has  a  thickness  about  twice  as  that  of  the  low  speed  layers. 
This  led  us  to  assume  a  bottom  model  comprising  of  three  layers  of  depths  8  m,  16  m  and  8  m 
respectively. 

In  the  analysis  of  data  in  the  frequency  range  250  Hz  -  3 10  Hz  there  was  a  problem  that  was  resolved 
by  re-assigning  the  mode  number.  Let  us  consider  the  model  M5.  The  bottom  model  is  a  two-layer  one 
with  top  layer  of  thickness  18  m  and  a  second  layer  of  thickness  12  m.  This  is  followed  by  a 
terminating  half  space.  The  compressional  wave  speed  in  the  top  layer  is  1700  m/s  and  in  the  bottom 
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layer  it  is  1550  /s.  In  the  terminating  half  space,  the  compressional  wave  speed  is  1850  m/s.  The 
density  is  1 .6  gm/cc  in  the  sediment  layers  and  in  the  half  space.  The  top  left  panel  of  Figure  4  shows 
the  spectrogram  obtained  by  performing  STFT  of  the  signal  at  the  receiver,  which  in  this  case  is 
located  at  a  distance  of  19270  m  from  the  source.  The  spectrogram  shows  that  the  mode  dispersion 
curves  for  all  the  modes  vary  smoothly  from  its  value  at  250  Hz  to  its  value  at  500  Hz.  Theoretical 
mode  dispersion  curves  for  Modes  11,  12,  and  13  obtained  using  the  known  water  column  and 
sediment  properties  are  overlaid  on  the  spectrogram  in  Figure  4.  The  mode  dispersion  curves  obtained 
from  theory  for  the  true  ocean  and  bottom  model  show  that  there  are  jumps  in  the  mode  dispersion 
curves.  Further  part  of  the  dispersion  curve  aligns  with  a  lower  mode  as  displayed  in  the  spectrogram. 
This  happens  for  resonant  modes.  It  is  seen  from  the  plots  in  the  top  left  panel  of  Figure  4  that  part  of 
Mode  13  as  per  theory  is  aligned  with  Mode  12  in  the  spectrogram.  Similar  features  occur  in  the  case 
of  other  modes  above  Mode  10.  The  incorrect  identification  of  modes  due  to  this  anomaly  will  lead  to 
large  errors  in  inverse  results.  In  the  top  right  panel  of  Figure  4,  we  plot  the  group  speed  of  Modes  1  to 
16  at  250  Hz.  These  are  represented  in  black  dots.  We  notice  that  two  values  (i.e.  those  for  Mode  1 1 
and  15)  are  much  higher  than  the  rest.  These  are  the  modes  where  the  resonance  is  observed.  We  now 
re-assign  the  mode  number  of  the  modes  eliminating  the  resonant  modes.  We  now  have  14  modes 
instead  of  16.  What  was  originally  Mode  12  has  now  become  Mode  1 1  and  so  on.  This  is  shown  as  red 
circle  in  the  figure.  We  do  this  for  modes  at  each  frequency  and  with  this;  a  modified  mode  dispersion 
data  is  obtained.  In  the  bottom  panel  of  Figure  4,  we  plot  the  modified  mode  dispersion  data  on  the 
spectrogram.  It  is  seen  that  the  curves  follow  the  modes  as  displayed  in  the  spectrogram.  In  adopting 
this  procedure,  we  identified  resonant  modes  i.e.  by  identifying  sharp  increases  in  the  group  speed  and 
eliminated  those  modes.  Simultaneously  we  made  suitable  changes  to  the  mode  function  and 
eigenvalues  associated  with  the  modes. 


Inversion  Result  Region  M2 


Figure  3:  Inversion  for  model  M2  without  discontinuity  incorporated  in  inversion  procedure  (red 
line).  The  black  line  is  the  true  model  and  the  blue  line  is  the  starting  model. 

The  results  of  the  inversion  performed  using  data  in  the  frequecy  range  50  Hz  -  1 10  Hz  and  in  the 
range  250Hz  -  3  lOHz  are  shown  in  Figures  5  for  four  of  the  seven  regions.  It  is  seen  from  the  results 
that  the  bottom  model  for  range  dependent  model  such  as  M2  that  are  closer  to  the  tme  bottom  model 
when  low  frequency  data  (i.e  data  in  the  frequency  range  50Hz  -  1 10  Hz)  are  used  when  compared  to 
the  results  using  data  from  the  higher  frequencies.  This  is  to  be  expected  as  mode  function  at  the  higher 
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frequecies  decay  fast  with  sediment  depth  and  hence  the  larger  errors  while  using  data  from  higher 
frequencies.  There  are  also  larger  errors  while  performing  range-dependent  inversions.  Bottom  models 
for  the  remaining  three  regions  also  show  similar  trends  in  the  results. 

E.  Discussion  of  Inversion  Results 

The  differences  between  the  true  bottom  models  and  the  models  obtained  by  inversion  of  the  mode 
dispersion  data  are  evident  in  Figure  5.  We  now  investigate  how  these  differences  in  bottom  model 
affect  acoustic  fields  predicted  by  these  sediment  models  even  though  the  differences  from  the  true 
model  are  small.  We  study  this  at  three  frequencies  50  Hz,  250  Hz,  and  500  Hz.  In  performing  this 
computation  we  considered  the  transmission  from  the  source  to  the  receiver  at  the  location  of 
Sonobuoy  2.  In  Figure  6  the  fields  as  predicted  by  the  inverted  bottom  models  are  compared  with  the 
field  as  predicted  by  the  true  bottom  model. 

The  field  at  50  Hz  as  predicted  by  model  from  50  Hz  -1 10  Hz  data  is  in  good  agreement  with  the  field 
predicted  by  the  true  bottom  model  up  to  about  5000  m  in  range  after  which  the  differences  increases 
with  range.  The  wave  number  spectrum  for  the  true  field  and  the  field  predicted  by  the  model  are  in 
very  close  agreement  except  in  the  case  of  Modes  3  and  4.  This  shows  how  even  small  differences  in 
the  wavenumber  spectrum  cause  differences  in  the  field  which  accumulate  with  range.  In  the  case  of 
field  obtained  using  the  bottom  model  from  250  Hz  -  3 10  Hz  data  the  fields  the  agreement  between  the 
fields  are  not  as  good  as  in  the  previous  case. 
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Figure  4:  The  top  left  panel  shows  the  spectrogram  with  the  theoretical  dispersion  curves  for  Modes 
11, 12,  and  13  overlaid.  In  the  top  right  panel  we  have  plotted  (black  spots)  the  group  speed  for  the 
first  16  modes  at  250  Hz.  Overlaid  are  the  group  speed  after  deleting  those  which  are  outliers  (red 
circle).  The  bottom  panel  shows  the  revised  dispersion  curves  obtained  after  deleting  the  outliers 

over  laid  on  the  spectrogram. 
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Compressional  Wave  Speed  (m/s) 


Compressional  Wave  Speed  (m/s) 


Figure  5:  True  bottom  model  and  the  inverted  model  for  bottom  models.  The  top  panels  shows  the 
result  for  Ml  and  M2.  The  bottom  panels  show  the  results  for  M3  and  M6.  In  the  plots  the  true 
bottom  mdoels  (Black)  are  compared  with  models  obtained  from  data  in  the  frequency  range  50Hz  - 
llOHz  (Red)  and  models  obtained  from  data  in  the  frequency  range  250Hz  -  310  Hz  (Magenta). 
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Figure  6:  (Row  l)The  plots  are  the  field  at  50  Hz  predicted  by  models  obtained  from  50  Hz  -110  Hz 
data  followed  by  field  predicted  using  model  from  250  Hz-  310  Hz  data.  (Row  2)  This  is  followed  by 
fields  at  250  Hz  predicted  using  models  from  50  Hz  -  llOHz  data  and  250  Hz  -310  Hz  data.  Beneath 
each  of  the  fields  are  the  plots  of  the  wave  number  spectrum  associated  with  these  fields.  (Row  3) 
The  fields  at  500  Hz  as  predicted  by  bottom  models  from  50  Hz  -110  Hzdata  and  250  Hz  -31 0  Hz 
data  are  shown.  Beneath  the  fieds  are  the  corresponding  wave  number  spectrum.  In  all  plots  the 
black  lines  are  the  predictions  by  the  true  bottom  model  and  the  dotted  red  lines  are  the  fields 

predicted  by  the  inverted  botom  model. 
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It  is  seen  form  the  wave  number  spectrum  that  Mode  4  in  very  weak.  Further  there  is  a  large  error  in 
the  wavenumber  of  Mode  3  between  that  of  the  true  field  and  the  field  from  the  inverted  model.  These 
differences  in  the  wavenumber  spectrum  contribute  to  the  differences  between  the  fields.  In  the  case  of 
predictions  at  250  Hz  the  predicted  field  generally  follows  the  true  field  in  terms  of  their  magnitudes 
though  there  are  differences  in  the  interference  pattern.  The  wave  nuber  spectrum  of  the  fields  show 
marked  differences  in  the  spectrum  which  results  in  the  differences  in  the  interference  patterns.  The 
same  is  true  in  the  case  of  field  predicted  at  500  Hz. 

CURRENT  WORK 

The  purpose  of  the  study  described  in  the  paper  is  to  investigate  whether  routine  Naval  operations  can 
be  suitably  modified  to  collect  data  that  can  be  used  to  estimate  sediment  properties.  Such  a  capability 
will  provide  critical  inputs  to  sonar  performance  evaluation  while  operating  in  shallow  water  areas. 

The  methodology  used  for  estimating  range-dependent  and  range  independent  sediment  properties 
have  been  well  documented  in  the  literature  [4,  5].  The  task  was  therefore  to  see  how  these  well  known 
method  will  perform  in  an  environment  which  has  complicated  range-dependent  sediment  structure 
and  at  different  frequency  ranges.  A  suitable  environment  was  created  using  the  well  surveyed  area  off 
the  coast  of  New  Jersey.  The  different  regions  of  this  environment  have  differing  depths  and  varying 
sediment  properties  were  assumed  for  these  regions.  The  inversions  were  carried  out  at  two  frequency 
ranges  (i.e.  50  Hz  -  1 10  Hz  and  250  Hz  -  310  Hz)  using  methodologies  developed  for  range 
independent  and  range-dependent  inversions.  The  results  of  inversions  show  that  it  is  possible  to 
extract  the  sediment  properties  though  the  errors  are  more  when  range-dependent  inversions  are 
involved.  The  evaluation  of  the  inverted  bottom  models  indicate  that  the  models  can  predict  fields  at 
50  Hz,  250  Hz,  and  500  Hz  but  with  errors  that  increase  with  range  as  expected.  Further  work  using 
data  from  higher  frequencies  namely  500  Hz,  1000  Hz,  and  above  is  in  progress. 

IMPACT/APPLICATIONS 

The  results  of  this  simulation  study  will  enable  validation  of  the  proposed  method  for  estimating  range- 
dependent  sediment  properties  using  data  collected  during  routine  naval  operations. 

RELATED  PROJECTS 

The  estimation  of  sediment  properties  in  a  shallow  water  environment  is  an  active  area  of  research. 

The  new  procedures  that  are  likely  to  come  out  of  this  work  will  therefore  be  of  great  benefit. 
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